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We use different particular classes of axially symmetric Szekeres Swiss-cheese models for the study 
of the apparent dimming of the supernovae of type la. We compare the results with those obtained 
in the corresponding Lemaitre-Tolman Swiss-cheese models. Although the quantitative picture is 
different the qualitative results are comparable, i.e, one cannot fully explain the dimming of the 
supernovae using small scale (~ 50 Mpc) inhomogeneities. To fit successfully the data we need 
structures of order of 500 Mpc size or larger. However, this result might be an artifact due to the 
use of axial light rays in axially symmetric models. Anyhow, this work is a first step in trying to 
use Szekeres Swiss-cheese models in cosmology and it will be followed by the study of more physical 
models with still less symmetry. 
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I. INTRODUCTION 

The Universe, as we observe it, is inhomogeneous al- 
most at all scales. However, the standard approach, 
based on the cosmological principle, assumes that homo- 
geneous models with Robertson- Walker (RW) geometry 
can be successfully employed in a cosmological frame- 
work, i.e., above some '-^lOO Mpc scale. However, within 
such a framework, one needs to assume that the Universe 
is filled with dark energy in order to achieve concordance 
with SNe la observations. But dark energy has never 
been directly observed neither in the Universe, nor in 
laboratories, and it has very exotic properties, namely, it 
is a kind of gas with negative pressure or a mere cosmo- 
logical constant with an amplitude too small to account 
for the vacuum energy in the standard model of particle 
physics. Therefore, it is of the utmost importance to test 
whether the observations can be explained with models 
that do not require dark energy but that instead take 
inhomogeneities into account. 

A number of approaches have been suggested so 
far. Among them the one patch spherically symmetric 
Lemaitre-Tolman models with a central observer [1-15], 
the thin shell approximations |16l4l9| , the Stephani mod- 
els [13-111, and the Szekeres models 0, [H. However, 
all these approaches have their restrictions. 

For example, for the Lemaitre-Tolman model assum- 
ing spatial spherical symmetry around the observer, this 
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symmetry has been turned into an argument stating that 
the observer's location at or near the symmetry center is 
too special as regards the requirement of the cosmological 
principle. Now, we want to stress here that these models 
should not be considered as assuming we are living at 
or near the center of any spherically symmetric universe. 
They should merely be considered as a first step in the 
process of modeling cosmological inhomogeneities, i.e., as 
a mere smoothing out of the inhomogeneities on angular 
scales, analogous to their smoothing out on the totality of 
space in homogeneous models. Thus, the subsisting inho- 
mogeneities are only radial. The use of such models must 
be merely regarded as a first approach which will be fol- 
lowed by more precise ways of dealing with the observed 
inhomogeneities. The Stephani class of solutions are con- 
formally flat [Ig and cannot therefore constitute general 
models for the Universe. The thin shell approximations 
are based on two RW models and are also spherically 
symmetric around the observer. The Szekeres models are 
the most general ones, even if the five functions defining 
a model depend only on the radial coordinate. 

An improvement is to use Swiss-Cheese models. Some 
models of this class where the holes were Lemaitre- 
Tolman patches have already been proposed in the lit- 
erature [7|. [27h31| . However, the results suggest that we 
need at least structures of several hundred Mpc in order 
to reproduce the supernova data. In Ref. [S^l the model 
proposed contains Schwarzschild regions, each with the 
mass of a galaxy. The results suggest that such a geom- 
etry cannot explain dark energy as an effect of inhomo- 
geneities. 

We propose here a generalization of these approaches, 
where instead of Lemaitre-Tolman or Schwarzschild 
patches, which are all spherically symmetric, we use Szek- 
eres models to represent the holes in the homogeneous 
cheese. The Szekeres model has several advantages, like 



anisotropic density distribution which leads to a different 
evolution pattern than in the Lemaitre-Tolman model 
|33L l34l |. Thus there is a possibility that the Szekeres 
Swiss-cheese models can come up with a more realistic 
picture, where the inhomogeneities are no more spherical 
and therefore can give a more physical representation of 
our Universe. However, for calculation simplicity and as 
a first step, we limit ourselves here to the study of an 
axially symmetric propagation of light through axially 
symmetric Szekeres holes. 

The structure of the present paper is as follows. Sec- 
tion|lI]is devoted to the presentation of the Szekeres mod- 
els and of the particular subclasses used in this paper. In 
Sec.lIIH we give the residual Hubble diagrams for the five 
axially symmetric quasi-spherical Szekeres Swiss-cheeses 
considered here and we compare them to the correspond- 
ing Lemaitre-Tolman models, to the ACDM model and 
to the actual data. In Sec. lIVi we display our conclusions. 



II. SZEKERES MODELS 

The Szekeres solutions 35] are the most general solu- 
tions of Einstein's equations one can obtain with a dust 
gravitational source. They have no symmetry, i.e., no 
Killing vector, and are therefore well-suited to describe 
a lumpy universe. Their metric in comoving coordinates 
and synchronous time gauge is 



ds' 



c^de 



^2a 



dr^ 



.2/3 



(dx^ 



dy% (2.1) 

to be deter- 



where a and /3 are functions of {t,x,y,r) 
mined by the field equations. 

There are two families of Szekeres solutions. The first 
family, where /3' = (here the prime denotes derivation 
as respect to r) is a simultaneous generalization of the 
Friedmann and Kantowski-Sachs models. Since it has 
found so far no application in cosmology, we do not dis- 
cuss it here. The second family of solutions is obtained 
when /3' ^ 0. When the Einstein equations are solved, 
its metric can be written, after a change of coordinates 
more convenient for our purpose [36[, 



ds' = c'dt^ - 



2 ($'-$i?'/i?)2 „ ^2(d:r2 + dy2) 



e-k 



dr^ - $^ 



E^ 



(2.2) 
where e — 0,±1, $ is a function of t and r and fc is a 



function of r. 



x-pV fy-Q''^ 



S 



s 



(2.3) 



with S{r), Pir), Q{r), functions of r. 



A. Quasi-spherical Szekeres models 

As can be seen from eq. (|2.2p . only e = -1-1 allows the 
solution to have all the three Friedmann limits (hyper- 



bolic, flat and spherical). This is induced by the require- 
ment of the Lorentzian signature of the metric. Since we 
are interested in the Friedmann limit of our model which 
we expect to become homogeneous at very large scales, 
i.e., that of the last-scattering, we will focus only on the 
e = -f 1 case. Such a case is called the quasi-spherical 
Szekeres model. 

Its metric, obtained for e — -1-1, becomes 



ds2 ^ c^dt^ 
where 

E=^- 



($' _ ^E'/Ef 



x-P 



$• 



(dx"^ + dj/2) 



£;2 



S 



(2.4) 
(2.5) 



In the spherical coordinates, E — S/{1 ~ cost?). 

Applying the Einstein equations to the metric p.4p 
and assuming the energy momentum tensor is that of 
dust, the Einstein equations reduce to the following two: 



1a2 2M , 1 _^2 
c2 $ 3 ' 



(2.6) 



where Af (r), is an arbitrary function related to the den- 
sity p via: 



Kpc 



2M' - 6ME'/E 
$2($' _$£;'/£;)■ 



(2.7) 



where k — 8nG/c^, A is the cosmological constant and 
Af is an arbitrary function of r. The 3D Ricci scalar is 



3 _ k ( ^k'/k-2^E'/E 

^$2 1^ $'-$S'/^ ^ 



(2.8) 



The Weyl curvature decomposed into its electric and 
magnetic part is 



E'^p = C".,psu^u 



7„* 



M(3$' - $M7M) , , 

— Vt ^^d ag 0, 2, -1, -1 , 



H, 



ap 



-gPa-^^iuC^" psu'u^ = 0. 



(2.9) 



As in the Lemaitre-Tolman model, the bang time func- 
tion, tB{r), follows from (|2.6p : 



d$ 



-k + 2A//$ + iA$2 



c[t-tB{r)]. (2.10) 



Since all the formulae given so far are covariant un- 
der coordinate transformations of the form r = g(r), this 
means that one of the functions k{r), S{r), P{r), Q{r), 
M(r) and tsir) can be fixed at our convenience by the 
choice of g. Hence, each Szekeres solution is fully deter- 
mined by only five functions of r: in the following, we 
choose fc, S, P, Q and M. 



B. Axially symmetric quasi-spherical Szekeres 
models 



As a first step, and for simplification purpose, we 
consider here axially symmetric quasi-spherical Szekeres 
models. Actually, since in our Swiss-cheese model de- 
scribed below, we use only radially directed light rays, it 
has been shown in [33| (see also [38J) that this implies 
that the Szekeres model should be axially symmetric. 
The simplest axially-symmetric Szekeres model obeys 



P{r) = xq ^ const, Q{f) 



yo 



const. 



In this case the dipole axis is along x = xq and y = yo 
(or in spherical coordinates along the directions "i? = 
and ■& = — tt). 

For the axially directed geodesies {dx = dy = 0), we 
obtain from (|2.4p 



dt _ 1 $' - <i>E'/E 
dr c y/l~ k 



(2.11) 



The plus sign is for r^ < Tq and the minus sign for r^ > ro, 
with Tg, the radial coordinate of the source and To, the 
radial coordinate of the observer. 

The redshift relation in this case is 13811: 



placed so that their boundaries touch wherever a light ray 
exits one inhomogeneous patch. Thus the ray immedi- 
ately enters another Szekeres inhomogeneity and does not 
spend any time in the Friedmann background. To match 
a Szekeres patch to a Friedmann background across a 
comoving spherical surface, r = constant, the conditions 
are: that the mass inside the junction surface in the Szek- 
eres patch is equal to the mass that would be inside that 
surface in the homogeneous background; that the spa- 
tial curvature at the junction surface is the same in both 
the Szekeres and Friedmann models, which implies that 
ksz = kpr^ and {ksz)' — "ikpr; finally that the bang 
time and also A must be continuous across the junction 
- unless otherwise stated we assume that A = 0. 

Besides matching the inhomogeneous patches we also 
need to take care of the null geodesies. However, as in 
this case we only consider axial geodesies the junction is 
trivial and requires only matching the radial, or equiva- 
lently the time component [39|. 

We consider 3 different types of models with 3 different 
backgrounds. The first model consists of small-scale in- 
homogeneities within an open FLRW model (i7m = 0.25). 
The others have larger patches within either an Einstein- 
de Sitter model or an open FLRW model with il^ = 0.9. 

Light propagation is obtained by solving eqs. (|2.13p . 
The luminosity distance is calculated as follows. First 
we solve for the angular diameter distance. Da, which is 
defined as [40!]: 



ln(l + z[ 
or equivalently: 



1 /",$'- ^E' E 
= ±- / dr ' , 

cj yr^fc 



dr c \/l — k 



dz 1 + z<^' -<PE'/E 

dt _ 1 $' - <PE'/E 
dz ^ 1 + 2;$' - ^E'/E' 



The Swiss Cheese models 



(2.12) 



(2.13) 



n2 



5S_ 



(2.14) 



where SS is the cross-sectional area of the bundle of null 
geodesies diverging from a radiation source, perpendic- 
ular to the propagation vector of light at a point with 
afEne parameter s and 5J7 is the solid angle subtending 
this area. The rate of change of SS is given by [4l[ : 



dSS 
"d7 



2ess, 



(2.15) 



where 6 is the expansion of the family of null geodesies. 
Using the Sachs propagation equations [4l|, |42l 



As seen above, in the case of axially directed geodesies 
the equations which describe light propagation simplify 
significantly. Moreover, density fluctuations (|2.7I) and 
curvature fluctuations [eqs. (|2.8p and (|2.9p ] are the 
largest along the axial axis {zLE' /E - when a light ray 
passes through the origin E' /E -> —E'/E, since E' /E is 
always multiplied by $ or <i> which are zero at the origin 
thus there is no discontinuity here.). Therefore in con- 
structing our Swiss-cheeses models we focus only on axial 
geodesies. 

When constructing a Swiss Cheese model, we need to 
satisfy the junction conditions for matching. Here we 
match the Szekeres inhomogeneity (holes) to the Fried- 
mann background (cheese). These Szekeres patches are 



dff 1 

^+e^ + \a\^ = - R^^k-k^, (2.16) 

ds 2 

^ + 20a = -ii?„^^,,e*"fc^e*A^r, (2.17) 

where Rap^v is the Riemann tensor, R^p is the Ricci 
tensor, |crp = (l/2)CTQ^cr"^, and a star denotes a com- 
plex conjugate. For axial geodesies e" = E{5'^2 + 
j^"3)/(v^$). Which implies that if initially cr = then it 
will remain zero. Therefore, the relation for the angular 
diameter distance is 



ds2 



Rank" kf^ Da. 



(2.18) 



Finally using the reciprocity theorem [40, |4j| the lu- 
minosity distance is 



Dl=Da{1 + zY 



III. RESULTS 



(2.19) 



A. Small scale (~ 50 Mpc) inhomogeneities 

We first study a Swiss-cheese which can be consid- 
ered as representing our local Universe with voids of 
size around 50 Mpc [4j|. Let us consider the Swiss- 



cheese model exhibiting holes each described by Szekeres 
patches whose 5 arbitrary functions, M , fc, Q, P, and S 
are of the following forms: 



M = Mo- 



Ml £3 
M2 exp 



{"^y 



-Ah{£-AXa)^ 



for e < Xa, 

for Xa ^ i ^ 3Xa 







for 3Xa ^ i ^ 4:Xa, 

ior e^4xa, 

(3.1) 

where I = r/kpc. Mo = (47rG/3c2)p,£3^ p, = "mS' 
a« = 0.25, Ho = 72 km s^^ Mpc-\ Xa = lO'^, Afi = 
a;-3M2e-l■^ M2 = -7 x 10" kpc. The above profile 
was chosen for the following reasons: it behaves like a 
FLRW model for £ ^ Xa but with lower density than 
outside, then for Xa ^ £ ^ 3xa we have a transition 
region, and a cubic behavior for 3xa ^ £ ^ 4:Xa, which 
allows for a smooth matching to the background values. 
Although the distribution of mass, i.e. the function M{r) 
is spherically symetric (it only depends on r) , the density 
distribution, as seen from l2.7l is not. For the same reason 
we chose the following profile for the function k: 

for r ^ Xh, 

for Xb ^ r ^ 3xb 

for 'Sxb ^ r ^ Axb, 
for r ^ Axb, 

(3.2) 

where fco = {upbC? /i — Hq)£'^, ki — k2X^^^e ^, fe — 

-8.84 X 10-6, xb = 8.68 x 10^. 




5=(5xl03 + £)±"-78 

P=l = Xo, 

Q = l = yo- 



(3.3) 
(3.4) 
(3.5) 



where -I- is for propagation from the origin [E' / E — 
0.78/(5 X 10^ + £)], and - towards the origin [E' / E = 
-0.78/(5 X 10^ -I- £)]. This model wih be referred to as 
model 1. As can be seen from (pT7)) - ([^ and (1511 - 
p.2p . for r > 40 Mpc model 1 becomes the homoge- 
neous FLRW model. The density profile along the axial 
geodesic is presented in the left panel of Fig. [T] First 
light propagates towards the center, E' > 0, and after 



passing through the origin, E' becomes negative, and so 
on. 

In model 1 we place the observer in the homogeneous 
region where Tq = 40 Mpc. We join Szekeres patches at 
r = 40 Mpc. Thus Szekeres patches are placed so that 
their boundaries touch and wherever a light ray exits 
one inhomogeneous patch it immediately enters another 
Szekeres inhomogeneity and does not spend any time in 
the FLRW background. 

The results for model 1 in the form of a residual Hubble 
diagram are presented in Fig. [2] The residual Hubble 
diagram presents Am as a function of redshift: 



Am 



TO — TO 



e7np 



5 log 



D, 



j-.e7np ■ 



where to*^"*^ and J^j^p are the expected magnitude and 
luminosity distance in an empty FLRW model. As seen 
in comparison with the corresponding Lemaitre-Tolman 
Swiss-cheese model the peaks are of larger amplitude. 
However, for z > 0.2 the fluctuations of the magnitude 
are negligible and the results are almost the same as in 
the hyperbolic FLRW background model. Thus as pre- 
viously noted [23, |30, |31| the Lemaitre-Tolman Swiss- 
cheese models with small scale inhomogeneities are not 
sufficient to explain away dark energy, and this is also the 
case for this particular axially symmetric quasi-spherical 
Szekeres Swiss-cheese. 



B. Medium scale (~ 500 Mpc) models with 
homogeneous center 

The Swiss-cheese model considered in this section pos- 
sesses holes that are described by the following functions: 



M = M< 



0- 



Mi£^ ior£i^Xa, 

Macxp -3(^)' for a:„ ^ ^ ^ 3a:, 

~Mi{£-4:Xaf iOT 3Xa i^ £ i^ 4:Xa, 

for ^ ^ Axa, 

(3.6) 

where £ - r/kpc, Mq = {AT:G/3c')pb£^ pb^^,Ho = 
72 km s-i Mpc-\ Xc = 1.65 x 10^ Mi = x-^M2e-^-^, 
M2 = -2.4 X 10^5 kpc, 



k ^ kn X 



ki£^ 
fc2 exp 



(^)^ 



kiiixd~£)^ 




for r ^ Xd, 

for Xd ^ r ^ 3xd 



for 3xd ^ r s^ ixd, 
for r ^ Axb, 

(3.7) 

where fco = {npbC^ /3 — H§)£'^, fci = k2X^^e ^, ^2 = 

-5.05 X 10-3 Xd = 1.46 X 10^ 



P=l = xo 



±el.5xl0~"r 



(3.8) 
(3.9) 
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FIG. 1: Density profiles along the axial geodesic for the models studied in this paper — model 1 (left), models 2 and 3 (center), 
models 4 and 5 (right). 
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FIG. 2: The residual Hubble diagram for model 1 (red). For comparison, the results for the corresponding Lemaitre-Tolman 
model (i.e. with M and k as in model 1 but with E' — 0) are also presented (green). The Union supernova data set [43| is 
depicted in the inset, for clarity only measurements for z > 0.2 are given. 



Q = 1 = yo- 



(3.10) 



where + is for propagation from the origin \E' jE = 
el. 5 X 10~^], and — towards the origin \E' jE — —el. 5 x 
10"*]. We consider two models: model 2 with e = 1 and 
model 3 with e = — 1. As can be seen from (|2.7|) - (j2.9p for 
r > 660 Mpc these models become the Einstein-de Sitter 
model. Also, from the construction it can be seen that 
for r < Xd — 146 Mpc these models are homogeneous. 

In models 2 and 3 we place the observer at the origin, 
where p/ Pcr ~ 0.28. We join Szekeres patches at r = 660 
Mpc. As in the previous case the Szekeres patches are 
placed so that their boundaries touch. Since we place 
the observer at the origin, we have two choices: we can 
send a light ray towards the largest dipole fluctuations 



{E' > 0) - model 2 
- model 3. 



or in an opposite direction [E' < 0) 



The results are presented in Fig. [31 As can be seen, 
for z > 0.6 the magnitude fluctuations have small am- 
plitudes. Moreover, the first peak is undesirable, since 
it omits all low-z supernovae, and, since low-z measure- 
ments have the lowest errors, the x^ is very high and is 
equal to 2419.79 for model 3 and 2291 for model 2. For 
307 supernova measurements this is an extremely high 



value^ . As in the previous case the resuhs are quah- 
tatively comparable with those reported when using the 
Lemaitre-Tolman Swiss-cheese models [2^, [2^ and show 
that these types of models are not suitable for fitting 
supernova data. 



C. Medium scale (~ 500 Mpc) models with 
inhomogeneous center and the minimal void scenario 

The Swiss-cheese model considered in this section ex- 
hibits holes that are described by the following functions: 



M 



PbKC^ 



df f 



Sp cxp - 



(3.11) 



where pb = ilm-^:^, O™ — 0.3, Hq = 68 km s"-^ Mpc~^, 
6p = 0.6 and a — 180 Mpc. The function k is calcu- 
lated from (I2.10p by assuming is = 0, and P ~ = Q, 
S = (10"^ -|- ^)^^°-*, where + is for propagation from the 
origin [E'/E = eO. 8/(10^ + £)] and - towards the ori- 
gin [E'/E = -e0.8/(103 + i)]. We consider two models: 
model 4 with s — I and model 5 with e = —1. As can be 
seen by using (P?7)) -(P ^ . for r > 400 Mpc these models 
become almost homogeneous. 

We place the observer at the origin. As in the preced- 
ing section we consider two cases: when light starts to 
propagate from the observer along E' > (model 4) or 
along E' < (model 5). We join inhomogeneous patches 
at r = 400Mpc. As seen by construction and from Fig. 
[1] (right panel) , this not a perfect matching — only the 
functions M and k are continuous along this boundary 
but not for example the curvature [eqs. (12. 8|) and ()2.9p ]. 
Thus, in addition to this not perfect matching, we study 
the minimal- void scenario, where there is only one in- 
homogeneous patch beyond which the model is almost 
homogeneous - see also [7| — since there is no matching 
here the results are exact. They are presented in Fig. 
m As in previous sections, the brightness fluctuations at 
hight redshifts are very small and thus for z > 0.3 the dif- 
ferences between the minimal void scenario (single path) 
and the Swiss-cheese models are small. Thus in terms 
of fitting the supernova data, there is no big difference 
whether the minimal void scenario or the Swiss-cheese 
model is considered. 

Model 4, because of the first peak which omits low- 
z supernovae, has very large x^, ie. x^ = 477 in the 
Swiss-cheese case and x^ — 517 in the minimal void sce- 
nario. The Lemaitre-Tolman model {E' = 0) presents a 
better fit, i.e. x^ = 277 for the Swiss-cheese model and 



^ X^ = Ei a I ™ where ^i and Ui correspond to the measure- 
ments of the 307 supernovae [45l | , ^o is the distance modulus in 
the empty FLRW model, and crint = 0.12. 



X^ = 279 in the minimal void scenario. The best results 
are obtained in model 5, where x^ = 269 in the Swiss- 
cheese case and x^ — 269 in the minimal void scenario. 
For the number of degrees of freedom of the sample ~ 307 
measurements - this is quite a good fit, but in comparison 
with the ACDM model less so. For this "concordance" 
model we get x^ = 249. However, if larger structures 
than ~ 500 Mpc are considered then one can obtain the 
same results as in the ACDM model (see for example [l^ 
where a model with a central Lemaitre-Tolman overden- 
sity is used to reproduce D^ with the same form on the 
observer's past light cone as in the ACDM model; or |l4| . 
the case of the cosmic fiow model, where the density at 
the current instant is homogeneous and the best fit model 
has x^ = 240). 



D. What if A / 0? 

So far we considered models with A = 0. In this section 
we would like to show that if A is included one can obtain 
better fits than in the standard case. 

Let us consider two models: one with small scale den- 
sity fluctuations; model 6 (which is basically model 1 
with flm = 0.25 but in addition we have JIa — 0.7) and 
one with large scale density fluctuations: model 7 (in 
this case a variation of model 5: but with Sp = 0.36, 
n,n = 0.25, riA = 0.6 and Hq = 72 km s'^ Mpc^^). 
The residual Hubble diagram is presented in Fig. [S] For 
model 6 we have x^ = 245 and for model 7, x^ = 239 (for 
the ACDM model x^ = 249). Not only these models have 
a better x^ but also the value of A changes compared to 
the standard one. In the case of model 6 the change is at 
the level of 5%, in the case of model 7 the cosmological 
constant changes by 15%. This shows that, in the era of 
precision cosmology, when we want to estimate the val- 
ues of the cosmological parameters with high accuracy, 
we cannot neglect the effect of inhomogeneities. 



IV. CONCLUSIONS 

We have used the quasi-spherical axially symmetric 
subclass of the Szekeres models to construct Swiss-cheese 
cosmological models which are generalizations of other 
Swiss-cheeses where the holes are spherically symmetric 
and that have already been proposed in the literature. 
This geometry has already proved to be very useful when 
studying structure formation as they allow for matter 
anisotropics |33ll34j. 

However, as regards supernova observations, we do 
not get qualitatively different pictures from the ones ob- 
tained with the corresponding Lemaitre-Tolman Swiss- 
cheese models. The results given in Sec . I IIII show that, 
as in Lemaitre-Tolman models [23, |30, [31] , small-scale 
inhomogeneities (~ 50 Mpc) do not alter the distance- 
redshift relation significantly, and thus cannot fully ex- 
plain away dark energy. To reproduce successfully the 
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FIG. 3: The residual Hubble diagram for model 2 (blue), model 3 (red) and for the corresponding Lemaitre-Tolman model 
(i.e. with M and k as in models 2-3 but with E' = 0). Supernovae measurements are extracted from the Union data set pSj . 
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FIG. 4: The residual Hubble diagram for model 4 (blue, SC - Swiss Cheese, MV - minimal void), model 5 (red, SC - Swiss 
Cheese, MV - minimal void) and for the corresponding Lemaitre-Tolman model (i.e. with M and k as in models 4 and 5 but 
with E^r ~ 0) are also presented (green). For clarity, the Union supernova data set [451 ] presented in the inset is for z > 0.2 
only. 



supernova dimniing, we must consider larger structures 
of order 500 Mpc, as in Lemaitre-Tolman Swiss-cheeses 
[7|. The presence of such large structures, in the frame- 
work of Lemaitre-Tolman models, can be ruled out by 
CMB constraints [49. However, for Szekeres models, the 



constraints might be different since, as shown in |39|, 
the CMB fluctuations are smaller than for the Lemaitre- 
Tolman Swiss-cheese models. 

The fact that in the studied cases we need large patches 
in order to explain supernova dimming without dark en- 
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FIG. 5: The residual Hubble diagram for models 6 (green), 7 (red), and ACDM (dashed line). For clarity the supernova data 
are presented for z > 0.1. 



ergy might be due to the axial symmetry of the models. 
Actually, it has been shown, that structure formation 
is five times faster in some non axially symmetric Szek- 
eres models than in the corresponding Lemaitre-Tolman 
models 33]. Since we know, from i28], that the evolu- 
tion of the inhomogeneities bends more the photon paths 



compared to the FLRW case than their mere geometry, 
it should be interesting to investigate if another geome- 
try of the holes might be able to enhance their evolution 
and therefore to have a stronger effect on the bending of 
the photon paths and, hence, on the supernova dimming. 
This will be the subject of future work. 
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